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ABSTRACT 

Even when we consider Newtonian stars, i.e., stars with surface gravitational redshift, 
z <C 1, it is well known that, theoretically, it is possible to have stars, supported 
against self-gravity, almost entirely by radiation pressure. However, such Newtonian 
stars must necessarily be supermassive(Hoyle and Fowler 1963; Fowler 1966; Weinberg 
1972). We point out that this requirement for excessive large M in Newtonian case, 
is a consequence of the occurrence of low z <C 1. But if we remove such restrictions, 
and allow for possible occurrence of highly general relativistic regime, z > 1, we show 
that, it is possible to have radiation pressure supported stars at arbitrary value of M. 
Since radiation pressure supported stars necessarily radiate at the Eddington limit, 
in Einstein gravity, they are never in strict hydrodynamical equilibrium. Further, it 
is believed that sufficiently massive or dense objects undergo continued gravitational 
collapse to the Black Hole stage characterized by z — oo. Thus, late stages of Black 
Hole formation, by definition, will have, z 3> 1, and hence would be examples of quasi- 
stable general relativistic radiation pressure supported stars (Mitra 2006) . It is shown 
that the observed duration of such Eddington limited radiation pressure dominates 
states is t ks 5 x 10 8 (1 + z) yr. Thus, t — > oo as Black Hole formation (z — > oo) 
would take place. Consequently, such radiation pressure dominated extreme general 
relativistic stars become Eternally Collapsing Objects (ECOs) and and the BH state is 
preceded by such an ECO phase. This result is also supported by our previous finding 
that trapped surfaces are not formed in gravitational collapse(Mitral 2005) and the 
value of the integration constant in the vacuum Schwarzschild solution is zero(Mitra2 
2005). Hence the supposed observed BHs are actually ECOs. 
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1 INTRODUCTION 

Stars may be defined as objects with an intrinsic self- 
luminosity which is generally sustained by the grip of self- 
gravity. Most of the stars that we know of are primarily 
supported by gas pressure p g rather than by radiation pres- 
sure p r . The usually insignificant role of p r is indicated by 
the parameter (Chandrasekhar 1967): 



where p — p g + p r is the total pressure. For most of the 
known luminous (Newtonian) stars /3 ~ 1. The actual role 
of p r is seen more clearly from the adjoint parameter: 

l-/3=£ (2) 
P 
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The relative role of p r may even be more directly expressed 
through the parameter 

_ Pr 1 - P 

x= — = —— (3) 

Pg P 

For real astrophysical objects, x would be highest at the 
centre and and in general decrease significantly towards the 
surface. The central value of x* = 0.006 for the Sun and 
the same for a 4Mq mainsequence star is 0.012 (as kindly 
pointed out by the anonymous referee) (0 denotes solar 
values). It may be borne in mind that the mean value of 
(1 — P) -C (1 — P*). It is known, however, that the value of 
(1 — p*) would rise significantly if the value of M would be 
very high though we hardly know of any single Newtonian 
star with M > 1OOA/0. In fact it may be pointed out that 
in any given stellar model, it is the value of (1 — /?*) which 
determines the mass of the star (Chandrasekhar 1967). Sup- 
pose one makes stellar models with a fixed composition (/x). 
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Then under the assumption of fixed chemical composition, a 
stellar mass of M will get defined from the equation (Chan- 
drasekhar 1967) (p. 82): 



1 



[3* 



M 

M© 



(4) 



so that M slowly increases with decreasing (3* (incresing x). 
Note M —* if (3* — > 1 or x* — ► in this non-quantum case, 
and to have a finite mass luminous star, one must have a 
finite value of x. Even in the case of (partially) degener- 
ate objects, a low value of (3, or a higher value of x would 
raise their masses. The usual Chandrasekhar mass, M c cor- 
responds to a perfectly degenerate object having ultimate 
relativistic degeneracy where momentum of the pressure giv- 
ing particles P — > oo and further the object is absolutely 
cold, T = 0. Moreover, it may be reminded that, the radius 
of Chandrasekhar 's critical white dwarf is zero, R c = 0! This 
shows that, only an infinitely dense singular configuration 
could be perfectly degenerate and have T — 0, and, on the 
other hand, all real finite configurations, in the strict sense, 
must be either partially degenerate or at a finite T. When 
the compact object has a finite temperature despite having 
assumed degeneracy, the mass upper limit gets modified as 
(Chandrasekhar 1967) (see p. 437): 



M c = M c (3 



-3/2 



(5) 



In order to keep the object almost "completely degenerate", 
Chandrasekhar imposed some artificial restrictions on the 
value of (3 which resulted in 



M c = 1.156M C 



(6) 



But the fact remains that perfect degeneracy, in a strict 
sense, corresponds to (3 = 1.0 or x = 0. It may be noted 
that the radius of this almost completely degenerate white 
dwarf too is zero, 1Z C = (Chandrasekhar 1967) (see p. 
441, eq. 140). However, theoretically, the possibility of hav- 
ing arbitrary finite x cannot be ruled out. As value of x 
would rise, the system would be more and more partially 
degenerate and soon, it would be more meaningful to call it 
"non-degenerate". In such a case, as (3 — > 0, one may have 
Mc — * oo, i.e., the very notion of an "upper mass limit" 
would vanish. Actually, the applicability of Eq.(5) would 
cease once degree of degeneracy would significantly reduce. 

As indicated by Eq.(4), in principle, an initially nonde- 
generate sequence of stars can always be turned into radia- 
tion pressure supported stars (RPSSs) by considering appro- 
priate higher value of M and vice-versa. Thus, irrespective 
of whether one follows the standard non-degenerate route 
or the degenerate route, one may, in principle, end up with 
non-degenerate radiation pressure supported stars. Now we 
focus attention on the question of mass range of such radia- 
tion pressure supported stars (RPSS) in the context of both 
Newtonian and Einstein gravity. 



2 GENERAL FORMALISM 

In reality, even if one would consider the stellar material 
to be pure hydrogen, at sufficiently high temperature, one 
would have pressure contribution from both electrons and 
protons. Also there will always be rest mass contribution 



by the electrons. But in the following, since we are inter- 
ested only in seeing some new qualitative result, we will 
consider an idealized fluid which is assumed to always re- 
main monoatomic. The total proper mass energy density of 
this fluid is 



Pa + (7 - 1) 1 Pg + Pr 



where 



pg = m p nc 



(7) 



(8) 



is the (baryonic) rest mass energy density, n is the proper 
baryonic number density, m p is proton rest mass, c is the 
speed of light, and 7 is the ratio of specific heats. Also 

p g = nkT (9) 

is the gas pressure, where k is Boltzmann constant and 

p r = aT A (10) 

is the radiation energy density. Here a = 7.56 x 10~ 15 erg 
cm~ 3 K ~ 4 is the radiation constant. As is known, 

Pr = \aT A (11) 

From Eqs. (9) and (11), we may also find that 

p g ink 

Further, instead of f3, it will be more convenient define a 
quantity f3 w = x _1 = p g /p r - One may also define a related 
parameter 

Pr aT 4 . . 

V = — = 2 ( 13 ) 

pg nipnc 1 

which is the ratio of the radiation energy density and the rest 
mass energy density. In Newtonian gravity, one has p r «/) Sl 
i.e., j<1, and, consequently, p » p g . In idealized models, 
fluids may be characterized by a polytropic equation of state 
(EOS): 



p = Kp 1+1/r ' 



(14) 



where K must be uniform within the fluid. If there would 
exist a strictly static self-gravitating configuration, where, 
K and m would be same everywhere, then one would have a 
polytrope of degree "m" . For the construction of stellar mod- 
els, one looks for polytropes which have a finite boundary 
radius R. Mathematically, this would mean that, p — p = 
for r > R. Both in Newtonian and GR case, in a strict sense, 
this vital boundary condition can be satisfied only if the object 
is absolutely cold T = so that p r — p r = for r > R and 
because one may indeed make p g = p g — at the boundary. 
In Newtonian gravity, since, p r = 0, and p — p g , this bound- 
ary condition can be satisfied much more accurately if not 
absolutely. In Newtonian gravity, one can consider the re- 
sults obtained by static polytropic models to be reasonably 
accurate. Now we shall attempt to qualitatively understand 
why Newtonian RPSSs must be of extremely high mass. 

2.1 RPSSs in Newtonian Gravity 

For the best known low mass Newtonian star, namely the 
Sun, as we know, the central value of f3, « f3 w « 0.003. 
Since T drops sharply as one moves away from the centre, 
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the mean value of fi must be extremely low in comparison 
to the above value. One may estimate the mean value of 
x for low mass Newtonian stars in the following way. By 
using Newtonian virial theorem for Newtonian stars, one can 
estimate the internal energy of the star as(Chandrasekhar 
1967) 



U 



GM 2 



(5 7 - 6) R 



(15) 



where 7 is the effective ratio of specific heats, G is the grav- 
itational constant and R is the radius of the star. As long as 
radiation pressure is really small, one can have 7 « 5/3. Fur- 
ther, if the mean temperature of the assumed monoatomic 
gas is T m , one has 

U=l—kT m (16) 
2 m v 



Then by combining Eqs. (15) and (16), it follows that 
_ 2 GMm v 



(17) 



7 R 

In terms of the surface redshift or compactness 2 = 
GM/Rc 2 , we can rewrite the mean temperature as 
2 

kT m = - z m p c (18) 

We may recall a that for Sun, 2 ~ 2 x 10 . This shows 
that the temperature of Newtonian stars are necessarily low 
because they are hardly compact, i.e., 2 <C 1. Alternatively, 
since, kT m is low in comparison to nucleon rest mass energy 
m p c , radiation pressure is low. In other words, radiation 
pressure is relatively low (for low M stars) because com- 
pactness is so low. One can ask then, what is the appropri- 
ate mean value of x for low mass Newtonian stars. To see 
this we invoke a relationship which shows that for z « 1, 
one has(Mitra 2006) 



y 

where 



GM 
a R^~ aZ 



a - 



L 

Led 



(19) 



(20) 



is the luminosity of the star in terms of its Eddington lumi- 
nosity L e d- By combining Eqs. (12), (13) and (19), we find 
that 



(21) 



3kT 



(22) 



Since we are interested only in the gross mean properties 
rather than precise central quantities, in the foregoing equa- 
tion, we may consider T = T m because virial temperature is 
a reasonably good measure of the gross mean temperature. 
Then we may use Eq.(18) in Eq.(21) to find that 

7a 

x « — x a 
6 

For Sun, L e d, at the surface is 1.3 x ~ 10 38 erg/s so that 
a ~ 3 x 10 -6 . Therefore the mean value of x for Sun is 
also ~ 3 x 10~ 5 . Again recall that, z m 2 x 10~ 6 for the 
Sun. Note the cosmic coincidence that, for Sun, the three 
apparently uncorrelated dimensionless quantities, namely, x, 
a and z have similar values within a factor of 10! In order 
to become a RPSS, the value of mean x has to rise from 



such an extreme low value to attain a stage of x 3> 1. Since 
for low mass stars, L ~ M /Br while Led ~ M, we will 
have, a = L/L e d ~ M 4,5 /i? ' 5 . For low mass main-sequence 
stars, crudely, M oc R, so that, a ~ M 4 . Thus if one would 
move to higher mass stars, initially, mean x would increase 
very rapidly. However, as x increases, Eq.(15) would need 
significant revision as the star would move away from being 
a gas pressure dominated 7^5/3 polytrope to a radiation 
pressure dominated 7^4/3 polytrope. In such a case, the 
pattern of variation of x with M would change drastically. 
It appears, that, in the regime of p r S> p g , the supposed 
Newtonian static RPSSs may be reasonably described by a 
Newtonian polytrope with m ~ 3(Weinberg 1972). Then, 
the mass of the star is(Weinberg 1972) 



M = 18M Q /3 U 



-2 2 



18M© x 2 fj 2 



(23) 



where fj, is the mean molecular weight. For a fully ionized 
hydrogen plasma = 2 It is interesting to invert Eq.(23) as 

^(^yv (24) 

which shows that in the large M range x increases relatively 
more rapidly as M 1 / 2 for the Newtonian stars. And eventu- 
ally, still, M has to be very high to ensure 1 > 1. To have 
a value of x ~ 10 (atleast) one must have, M ~ 1800Mq /i 2 
which is already very high. Note that, since, in principle, x 
can be arbitrarily high, again, there is no upper limit on the 
mass of self-gravitating configurations once we allow exis- 
tence of sufficiently high radiation pressure. We may also see 
how the "compactness" of the Newtonian star would change 
in this very high M range. To appreciate this point further, 
we consider specific model of Newtonian RPSSs( Weinberg 
1972) (see Eq. [11.5.9]): 

R oc /3~ 2/3 oc x 2/3 (25) 



Using Eqs. (23) and (25), we find then, that, in Newtonian 
case 

z oc M/R oc x 4/a (26) 

Further, using Eq.(24) into the foregoing equation, we see 
that 



z oc (M 1//2 ) 4 ^ 3 oc M 2 ^ 3 



(27) 



also increases quite rapidly with mass. Thus compared to 
the low mass Newtonian stars, in the high mass range, on a 
relative scale, one can have quite high values of z. However, 
this does not mean that for Newtonian stars one can ever 
have z 3> 1. This point has been emphasized by Weinberg: 
• The entire derivation leading to the concept of any 
Newtonian star necessarily assumes y = p r /p g <C 1. And 
this assumption must be ensured for (Newtonian) super- 
massive stars as well. As shown by Weinberg, in such a case, 
one must always have(Weinberg 1972) 

z = Tjjjs" « °' 39 ( 28 ) 

In summary, the fundamental reason that Newtonian RPSSs 
stars are extremely massive is that, by definition, and, 
in particular, z <C 0.39. It may be also pointed out that 
since in general L oc M d , where, d > 1 (for low mass stars, 
d ~ 5), and L e d oc M, a in general increases with M. Then 
Eq.(27) would show that, atleast, in the Newtonian regime, 
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in general, a increases with z. However, the maximum value 
of a — 1, and all the RPSSs, Newtonian or Relativistic, have 
this maximal value of a because a larger a would disrupt 
the star due to excessive radiation pressure. 



3 RADIATION PRESSURE SUPPORTED 
STARS IN EINSTEIN GRAVITY 

The general definition of "compactness" in GR may be given 
in terms of the surface redshift 



z = (1 - 2GM/Rc 2 )~ 1/2 - 1 



(29) 



One then easily finds that when GM/Rc 2 <C 1, (i.e., in the 
truly Newtonian regime), one has 



GM/Rc 2 



(30) 



It also follows that, the Event Horizon of a Black Hole (BH), 



defined by R 



Rs 



2GM/c corresponds to z 



And this explains why gravitation is extremely strong for 
BHs. As is well known, very massive objects undergo con- 
tinued collapse to become Black Holes. Thus, by definition, 
very massive objects, during their continued collapse, must 
pass through stages having z > any finite number. Al- 
though, it is not necessary, we may nonetheless mention that 
the exterior spacetime of any contracting self-gravitating 
object is represented by radiating Vaidya metric(Vaidya 
f951) which allows the possibility that z — > co. Further, it 
has been shown that, whenever, self-luminous objects have 
z > 1 (Mitra 2006), one will have 



y ~ a z/2 



(31) 



Then by combining Eqs. (12) and (13), Eq.(22) gets modi- 
fied, in the extreme GR case, as 



a (z/2) 



3kT 



(32) 



In all physically realistic cases of self-luminous objects, a 
is always finite. In the Newtonian case of z -C 1, we 
found that, a increases rapidly with M and hence with 
z. Since increase of z implies stronger self- gravity, we may 
say that as if stronger gravity churns out more radiation 
from self-luminous self-gravitating objects. And as z — > oo, 
the entire object becomes a ball of radiation/pure energy 
(pr/pg — ► oo). Note that this happens before the forma- 
tion of any Event Horizon(EH) and which indicates that 
the EH is actually synonymous with the central singularity 
and the integration constant of the vacuum Schwarzschild 
solution has the unique value of zero(Mitral 2005; Mitra2 
2005). Actually, one can have continued gravitation collapse 
for arbitrary (small) M too provided p would be suitably 
high. In contrast to the Newtonian case, where a (relatively) 
higher z demands, higher M (Eq.[27]), in GR, the very fact 
that there could be continued gravitational collapse at any 
mass scale means that one can have high z >• 1 for arbi- 
trary mass and hence, one may have a — ♦ 1, at suitably 
high value of z, at arbitrary mass scale. For instance, for 
the Sun, a « 3x ~ 10~ 5 and z ^ 2 x 10~ 6 . But if we con- 
sider a collapsing hot proto-neutron star in its final stage 
when it is giving birth to a hot neutron star, a w 10~ 3 as 
z ~ 10 _1 (Mitra 2006). Note the rise in the value of a with 



z even at a low mass scale of M ~ IMq as we are about to 
enter the regime of Einstein gravity. 

The gravitational mass energy of the star is defined as 



Mc 



-f 

Jo 



4nr p dr 



The mean value of mass energy density is thus 
3Mc 2 



4ttR 3 



(33) 



(34) 



From Eq.(29), note that in the 2>1 range, the radius of 
the contracting body would be hovering around its instanta- 
neous Schwarzschild radius, i.e., R « R s = 2GM/c 2 . Using 
this fact in the foregoing equation, we have 



P = 



3c 8 



32ttG 3 M 2 



(35) 



Note that Eq.(31) implies that, for z > 1, p r > /) s so that 
total p = p r + p g « p r . Further recalling that p r — aT 4 , 
from Eq.(35), we obtain 



T = 



3c s 



32tt a G 3 M 2 



1/4 



(36) 



Numerically, one finds, that, for such a state, the mean tem- 
perature of the body in this phase is : 



600 



M 



-1/2 



MeV 



By substituting Eq.(37) in Eq.(32), we obtain 
M ' 



15M<7 



(37) 



(38) 



Note that, even in the extreme GR case, x increases with 
M, and, in particular the behaviour x oc M 1 ^ 2 is just what 
we found for Newtonian RPSSs. But unlike the Newtonian 
case, in the GR case, one can have z > 1. Eq.(38) shows 
that, with a > 0, and, z — > oo, one would be able to satisfy 
the condition x 3> 1 at arbitrary mass scale for appropri- 
ate high value of z 3> 1. As mentioned earlier, occurrence 
of x 3> 1 means a — 1 (rather than a S> 1, which would 
disrupt the star). For isolated bodies, radiation pressure and 
energy density are directly associated with outward radia- 
tion/heat flux. This outward radiation flow has a repulsive 
action on the plasma. As radiation pressure tends to grow 
unabated with unabated increase of z, repulsive effect on the 
plasma must be able to counterbalance the pull of gravity 
at some stage. And this stage corresponds to attainment of 
L ~ L e d or a m 1. In fact, domination of radiation pressure 
is a much less demanding phenomenon than domination of 
radiation energy because recall that though the Newtonian 
supermassive stars have p r 3> p g , they still have p r <C p g . 
On the other hand, since Eq.(31) demands p r 3> p g at ap- 
propriately high z, such an occurrence actually automati- 
cally denotes domination of radiation pressure. Therefore, 
for x S> 1, we may rewrite Eq.(38) as 



M 



\15Mq 



1/2 



0.25 z 



M 



1/2 



Z > 1 



(39) 



The absence of z in Eq.(24) demands that one can have 
i»l only for very large values of M. On the other hand, 
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the presence of z in Eq.(39) shows that, in the extreme rela- 
tivistic case, one can have RPSSs for an arbitrary stellar 
mass, high or low. Hence as one would proceed towards 
z — > oo during continued gravitational collapse, one must 
obtain radiation pressure supported configurations (x ^> 1) 
at arbitrary mass rather than finite mass BHs. Although we 
worked only in the regions of z <C 1 (Eq.[24]) and z » 1 
(Eq.[39]), the close similarity in the forms of Eqs.(24) and 
(39) strongly suggests that Eq.(39) may be valid for RPSSs 
in the entire relativistic range range of z > 1. 



4 DISCUSSION 

Newtonian stars are defined by z < 1 as well as p r <C p g , 
i.e., y -C 1. It is such self-imposed constraints which cause 
x <C 1. However, even in the Newtonian regime, i.e., despite 
having 1/ < 1 and z -C 1, one can have RPSSs (x 3> 1) 
for masses above 7200Mq. Hence we find that the non- 
occurrence of low mass RPSSs in the Newtonian regime is 
intricately linked with the occurrence of z < 1 and y <C 1 
in such cases. On the other hand, GR is unleashed in full 
glory during continued gravitational collapse when one can 
have z > 1 and 1/ > 1 for an arbitrary mass, low or high. It 
naturally followed then that, in the extreme GR case, one 
can have a RPSS even at arbitrary low M. While occurrence 
of Newtonian supermassive stars may be only a theoretical 
possibility, Einstein RPSSs must be a reality because of the 
following simple reason: 

As z tends to increase indefinitely during continued 
collapse, the strong gravity almost completely traps the 
collapse generated neutrinos and photons within the body 
of the star(Mitra and Glendenning 2006). The density of 
trapped radiation also increases because of stellar mat- 
ter -radiation interaction, i.e., the diffusion of the internal 
radiation(Mitra 2006). As a result, the trapped radiation 
pressure and energy density increases at least as fast as 



jca _ 



R^(l 



-(Mitra and Glendenning 2006). On 



the other hand, the locally defined Eddington luminosity 
grows as(Mitra 1998) 

4itGMc , 



Led — 



(40) 



where k w 0.4 cm 2 g 1 is the Thomson opacity. The radia- 
tion pressure associated with L e d is 



1 L e 



Ped = - - 



oc R (l + z) 



(41) 



3 4tt.R 2 c 

Therefore, p r /p e d oc i? _1 (l + z). Thus as R decreases and 
l + z increases p r — > p e£ j eventhough, initially, p r <C p e d and 
L <C L e d. It is at this stage that L — > L e d and a — > 1. Si- 
multaneously, one will have both p r ^> p g and p r ^> p g at 
this stage. Specific models of GR continued collapse which 
do not restrict p r by any means do clearly show that repul- 
sive effects of unabatedly rising of p r and p r may not only 
counterbalance the inward pull of gravity but even cause 
the fluid to bounce back(Herrera and Santos 2004; Herrera, 
Prisco & Barreto 2006). Such studies are fully consistent 
with the generic picture of formation quasistatic GR RPSSs 
at arbitrary mass scale discussed in this paper. The observed 
luminosity of such RPSSs will however be lower by a factor 
of (l + z) 2 than the local value because of joint effect of 
gravitational redshift and gravitational time dilation: 



A-kGMc 



(l + z) 2 k(1 + z) 



(42) 



Consequently, even if the system would be assumed to be at 
a given fixed value of z = z, its observed duration as seen 
by a far away astronomer would be 



t(z) 



Mc 2 



kc 
4^G 



(1 + *) 



(43) 



Since in principle, during continued collapse, z — > oo, clearly, 
the observed time scale for depletion of mass energy be- 
comes infinite for arbitrary value of the opacity k: t = oo. 
Hence the RPSSs tend to collapse for infinite duration in 
order to attain the BH (z = oo) state and, therefore, may 
be called as "Eternally Collapsing Objects" (ECOs) (Mitra 
2000, 2002). It has been also shown that since the eventual 
BH mass would be zero, the comoving proper time for its 
formation would also be infinite(Mitra 2002). Since the ob- 
served BH Candidates must be formed in gravitational col- 
lapse and of finite age, they must be ECOs (z 3> 1) rather 
than true BHs (z = oo). In retrospect, long back the RPSSs 
were suggested as the central engine of quasars(Hoyle and 
Fowler 1963; Fowler 1966). However this attempt failed be- 
cause such RPSSs (i) were considered to be either Newtonian 
or Post Newtonian objects with low temperatures, (ii) The 
basic source of energy liberation was considered to be of nu- 
clear origin. In contrast the relativistic RPSSs considered 
here are fed by energy release due to secular gravitational 
contraction and the source of energy is the entire mass en- 
ergy (E — Mc 2 ). Even if they would momentarily be un- 
stable, the contraction generated luminosity would ensure 
that they pass from one quasistatic state of z = zi 3> 1 to 
another with 22 > z\. 
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